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1. Computation of the M Factor 


Faraday rotation measurements between station and satellite are 
affected by both the earth’s magnetic field and the ionosphere, but can 
be reduced with the aid of proper conversion factors to a measure of the 
ionosphere alone. The INTASAT satellite transmits plane- polarized 
signals at 40. 01000 and 40. 01025 MHz, These frequencies are much 
higher than the electron collision frequency and the gyro- and plasma 
frequencies in the ionosphere; thus, a ’quasi- longitudinal’ approximation 
will hold for propagation in all directions making angles of less than 
about 89- 5° .with the earth's magnetic field. Using a simplified form of 
the Appleton-Hartree formula for the phase refractive index, a relation- 
ship can be obtained between the Faraday rotation angle along the angular 
path and the total electron content along the vertical path) intersecting 


the angular at the height of maximum electron density. 



cos 0 Nds = 


K 

f 2 


f hu 

\ B cos 0 sec y Ndh 


( 1 ) 


0 = Faraday rotation angle in degrees 
K = 1. 699 ~ constant 
f = signal frequency in hertz 

B = earth's magnetic field strength in ampere-turns/m 
0 '= angle between direction of propagation and earth's magnetic field 
y = zenith angle 

N = electron density in electrons /m 3 
s = path length in m 

h = height above surface of earth in m 

h u = upper integration limit is the height of the INTASAT satellite 


Using the second mean value theorem of integration, the function 

B cos 0 sec y is removed from under the integral sign and replaced 

by a ’mean* value. 

v r h u 

n = -pr M V Ndh = M N t 


0 
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( 2 ) 



M = ‘mean’ value of (B cos 8 sec x) i n ampere-turns/m 
N t = vertical total electron content in electrons /m 2 column 


The conversion factor M is obtained from both of the above expressions 
for 0, as, r h, 


'• as > rh u 
M = Q - — 


B cos 0 secX'Ndh 


j*h u 


( 3 ) 


Ndh 


The integrals are evaluated in computer mode by generating the electron 
density N and the function (B cos 0 sec X N ) at various height intervals 
and numerically integrating by Gaussian quadrature. The eLectron density 
at each height h. is calculated by the worldwide Bent Ionospheric profile 
model (Refe rence 1 &:2), Each parabolic and exponential segment of the profile 
was integrated separately with a varying number of points to achieve 
maximum accuracy. A total of 23 points was used to evaluate the integrals 
defined in equation (3). The components of the magnetic field strength 
are obtained by a spherical harmonic analysis routine as described iri 
Appendix B. The assumption of straight line propagation through a 
spherica ily stratified ionosphere was made. No bending corrections were, 
calculated as this would have required a prohibitive amount of computer 
time. At the INTASAT fr equencie s, bending is a second order effect. 

Given the straight line propagation as sumption the zenith angle at each 
height h then becomes a function of the ground elevation angle, and the 
angle 0 is calculated using the station and satellite positions and the 
direction of the magnetic field. 


2 . Computer Listing of the M Factor 

The M factors are printed on the computer listing for 3 9 station receiv- 
ing signals from the INTASAT satellite during the specified time period. 

The data is sorted by station and date. 

For each day the visible satellite passes are numbered sequentially 
starting at one. If the satellite is continuously visible past 24 hours, the 
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last pass of the first day will only be partial. However, the first pass of 
the following day will list the complete pass, repeating the data from the 
first day and flagging the time column by * to indicate the day change. 

The Greenwich Mean Time for each day runs from 0 hours 0 minutes 
0 seconds to 23 hours 59 minutes 59 seconds. Time values of 23:59:59. 5 
or greater, but less than 24:00:00 are rounded to 24;00:00. 

The ionospheric pierce point is printed as the latitude and longitude 
at which the angular ray passes through the maximum electron density 
along the path. At this location, the ionospheric profile is computed by 
the Bent model as required for the computation of M. The M factors are 
listed in units of ampere-turns/m, and related to Gauss units by 1 Gauss = 

79. 58 ampere-turns /m. If the M value is flagged by the angle 0 between 
the direction of propagation and the magnetic field has obtained values 
between 89, 5 ^0^ 90.5°, for which the equation relating the 1 Faraday 
rotation and the total electron content is no. longer valid. If this condition 
occurs above 1000 km, an estimate for M is computed using the same 
equation; if the condition occurs below 1000 km, however, M is not 
computed and a zero value is printed. 

Total vertical electron content N T (el/m 2 } is reduced from the Phraday 
rotation measurement 0 (deg. ) using the M factor (amp-turns /m) by, 

KD, 

N t = f 2 M , . (4) 

where f is the signal frequency (Hz) and K=l. 699 is a constant. 

An example of the computer listing is given in Appendix C. 

3. Variation of the Faraday Factor 

A number of graphs are included to demonstrate the variation of the 
Faraday factor with local time and season, with magnetic latitude, elevation 
and azimuth angles. The effect of typical day to day fluctuations on the 
Faraday factor due to sudden increase and decrease in the ionospheric 
density and height are shown as well as the changes in the angle between 
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the direction of propagation and the magnetic field lines. 


As frequently used for convenience, the Faraday factor F in 
the Figures is the quantity computed from, 

N t (e/m 2 ) = FO{degrees ), (5) 

giving the direct conversion for the angular measurement to the vertical 
content for a signal frequency f=137 MHz. The relationship to M is given 


b Y» 


M (amp. -turns /m) = 


J 2 _ 1. 105X10 1 6 

KF " F(1 /m 3 degrees) 


( 6 ) 


Figures 1 through 5 point out the importance for modeling the Faraday 
factors correctly with respect to the station position, where the magnetic 
latitude is of most significance, and with respect to the direction of obser- 
vation, since the elevation ‘and azimuth angles determine the direction at 
which the magnetic field lines are intersected as well as the location at 
which the wave passes through the densest part of the ionosphere. Less 
important are the specific season and diurnal influences producing variations 
of only about 2 to 6% in the Faraday factors, as well as the day to day 
prediction errors in f c F2, having even less effect. However, prediction 
errors in ionospheric height which could easily be caused by sudden day 
to day changes, can have a significant influence on the Faraday factors 
especially for obse rvations along angular paths. Variations of 4- 100 km 
in height are not uncommon particularly in the equatorial region. Errors 
of 5% in the Faraday factor are typical for paths at vertical incidence, 
but as shown in Figure 3b. for angular paths errors of around 3 0% in 
the Faraday factor might occur resulting in proportionally large errors 
in N t , since N T =FO. The predicted values of the height of maximum 
electron density obtained from the Bent model are on average within the 
accuracy of the measured values, which considering instrumental and 
reduction techniques, are about 15 km. However, the day to day variations 
are quite a bit larger, and on occasion, deviations in the predictions of 
100 km from the height measurements have been noted. 


- 4 - 


For a number of stations and observation angles Figures 6a- e. demon- 
strate the behavior of the angle 6 between the direction of propagation 
and the earth's magnetic field lines between heights of 100 and 1000 km. 

For fixed station positions' and elevation angles the 0 angle versus height 
curves are shown for various azimuth directions. When the condition 
89. 5 0 ^ 90. 5 0 occurs, the equation relating the Faraday rotation angle 

and vertical electron content no longer holds true. When 0 passes through 
90° at a certain height, the wave experiences rotation of the polarization 
vector in one direction from the satellite down to that height, and rotation 
in the opposite direction below that height. Contributions to the rotation 
of the polarization vector in reversed directions cancel out, thus the , 
measurement is not representative of the ionosphere between the satellite 
and the station. 
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Figure 1. Seasonal and Diurnal Variatio 
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Figure 5a. Variation of the Faraday Factor with Changes in Elevation 
and Azimuth Angles at 80° Magnetic Latitude. 

Station Position - 68.6°, 279.4% Date = 16 March 1967, UT = 12 hour 
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Figure 5b. Variation of the Faraday Factor with Changes in Elevation 
and Azimuth Angles at 3 9° Magnetic Latitude. 

Station Position=28. 6°, 279. 4°, Date = l6 March 1967, UT=11 hours 
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Station Position =-1.2°, 279.4°, Date = 16 Mar 1967, UT = 14 hours 
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APPENDIX A 


DESCRIPTION OF THE BENT IONOSPHERIC MODEL 

A: 1 I onospheric Model Development 

For several years scientists have investigated many different 
approaches to modeling the ionospheric profile on a theoretical basis. The 
names and types of these methods are well known and will not be discussed here, 
but it is obvious after all the years that a good theoretical ionospheric profile 
still does not exist. 

The object of our past investigations was to come up with an ionospheric 
profile that could give much improved results for refraction corrections in 
satellite communications to ground or to another satellite than had been obtained 
with the Chapman and many other theoretical profiles. It would have been 
pointless for us to sit down and investigate another theoretical approach when 
so many more competent scientists are working on this problem. For this 
reason we decided that in this present time of computers, an empirical model 
taken from a vast data base may provide us with the profile we were looking for. 

It was our intention to acquire ionospheric data of any kind that helped 
us build up a data base covering minimum to maximum of a solar cycle and 
providing information up to 1000km. The lower layers of the ionosphere were 
neglected in terms of their irregularities although their electron content was 
added into the larger F layer; this was done to simplify the approach and as 
the prime objective was to obtain refraction corrections through the ionosphere, 
or at least to a point above 150 km, such an elimination would not be very 
detrimental. r 

Data from bottomside ionospheric sounders was obtained over the 
year 1962 through 1969 covering 14 stations approximately along the American 
longitudes having geographic latitudes 76 degrees to -12 degrees or magnetic 
latitudes 85 degrees to 0 degrees. This data was in the form of hourly profiles 
of the ionosphere up to the f Q F2 peak. Topside soundings were acquired for 
the years 1962 to 1966 covering the magnetic latitude range 85 degrees to 
-75 degrees and providing electron density profiles from about 1, 000 km down 
to a height just above maximum electron density. As the topside data was 
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®not available near the solar maximum, electron density probe data was 
obtained from the Ariel 3 satellite over the period May 1967 to April 1968 
from 70 degrees north to 70 degrees south geographic latitude and linked 
in real time to foF2 values obtained from 13 stations on the ground. 

A.l. I Ionospheric Profile 

In order to analyze the vast amount of data that was obtained a 'number 
of assumptions had to be made. . In the first case the topside sounding data 
did not geographically cover the entire globe and the bottomside data was 
only available for land masses and not over the oceans; however, as a local 
time effect is far more significant than a longitude effect, the data was 
analyzed as a function of latitude and local time. Geographic longitude was, 
however, taken into account for the determination of maximum electron density 
by using the ITS coefficients for f 0 F2 which are a function of latitude, longitude, 
time and solar activity. Secondly a theoretical profile was determined to which 
the data would fit. This profile which is used in the evaluation discussed later, 
is shown in Figure 7 and is the result of earlier work by Kazantsev (Reference 4), 
and unpublished work of Bent (1967) while at the Radio and Space Research 
Station in England and requires the knowledge of the parameters k : >k 2 ,k 3 , 

Yt , Y » f 0 F2,. and h B . The equation of the upper topside is exponential, namely, 

N = ,N 0 e" ka , 

the lower ionospheieis a bi-parabola, 

N = N. (l . iL) 2 , 
y H 

and the top and bottomside are fit together with a parabola, 

N = N - 0 - £.) • 

y? 
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where. 


N is the electron density 

N o is the maximum value of electron density 

N Q is the maximum electron density for each exponential 

layer 

a and b are vertical distances 

y^ is the half thickness of the lower layer 

y t is the half thickness of the upper parabolic layer 

k is the decay constant for an exponential profile.;' 

The upper parabola extends from the height of the maximum electron 1 
density up to the point where the slope of the parabola matches the slope of 
the exponential layer. The data investigated included over 50, 000 topside 
soundings, 6, 000 satellite electron density and related f 0 F2 measurements, 
and over 400, 000 bottomside soundings. 

A, 1,2 Topside Ionosphere 

The initial approach was to take the topside soundings and break them 
down into zones 5 degrees of latitude by 40 minutes of local time eliminating 
data in the same zones that have similar times and profiles, and therefore 
are duplicated. This resulted in over 1, 200 different areas in the northern 
and southern hemisphere with a reasonably constant density of data in each 
area. By these means it was possible to investigate the decay constant k 
in the exponential topside profile as a function of local time, latitude, solar 
flux, sunspot number and season. One of the major concerns was whether 
the decay constant k would be uniform for each sounding over the range 
1, 000 km to the minimum height, and investigations showed that such an 
exponential profile does not exist. The layer was, therefore, divided into 
three equal height sections from 1, 000 km to the minimum recorded height 
and the exponent k computed for the center point in each section. Figure 7 
shows such a division where the values under investigation are the decay con- 
stants kj, kg, kg. In most cases the topside soundings do not reach the height 
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of maximum electron density and therefore the gradient at this lower point was 
mathematically equated to the point where the gradient of the 'nose' parabola 
was the same. Extensive analysis of the acquired data showed these gradients 
to be similar, on average, at a height y. /4 above the maximum. electron density. 
At this point the value of f k F2. which defines the lowest point of the topside 
sounding, is 0. 93 f 0 F2. (N 0 in Figure T is the equivalent electron density to 

the frequency f k F2). 

For ah initial test the decay constants k for each of the three layers, upper, 
middle, and lower topside were plotted as a function of magnetic latitude and 
f k F2. Values from the northern and southern hemispheres were treated indepen- 
dently at first, but the analysis showed that there was excellent correlation 
between the two. Figure 8 shows the relationship between the three decay 
constants k and magnetic latitude for alllocal times, solar activity, and season. 
The equatorial anomaly and a 40 degree trough show in the lower topside layer. 
The 65 degree trough is not as evident as it is when the same analysis is done 
for various local times which suggests the physical variances of these anomalies 
should be investigated in more detail. 


It was found that correlations in k for specific f k F2 did not bear any 
further local time correlation, but bore a significant variation with solar 
activity and magnetic latitude. However, the correlation with solar flux was 
considerably better than that with sunspot number, even allowing for the delay 
in the effect reaching the ionosphere, so all further correlations were with 
the Ottowa 10. 7 cm solar flux. All these correlations were then plotted in 
graphical form to enable firal interpolation. 

Unfortunately the Alouette data did not cover the period at the peak of 
the solar cycle, but the Director of the U.K. Radio h Space Research Station 
made available electron density data from the Ariel 3 satellite to cover this 
period. The data had already been reduced thoroughly and the satellite electron 
density at about 55 0 km was provided with the sub- satellite f Q F2 value obtained 
from 13 stations around the world. If the satellite was not directly over an 
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ionosonde at the time of observation, the fc>F2 values from two or three 
transmitters in the general area had been interpolated in time and position to 
give the sub- satellite value. These interpolations had been carried out taking 
care to modify the values for uneven ionospheric gradients. Data that was in 
doubt was eliminated. While these values did not give the three exponential 
decay constants at each point, it was found that for similar conditions of solar 
flux and position, the Ariel 3 data fit very closely to the profiles deduced 
from Alouette 1. The profile equations developed for the lower solar activity 
period related to the topside sounders could, therefore, be extended to the 
larger solar flux values and still he in good agreement with the Ariel 3 data. 
Typical results from this analysis are shown in the graphs of Figure 9. The 
original data curves were less regular, and since the variations were mainly 
caused by the relatively low data density in each group after division of the 
large data base, the data was smoothed by the fitting of straight lines. In 
order to interpret these graphs and obtain a profile, we need the value of f Q F2. 
and the magnetic latitude position. These values will indicate which gro-ph 
relates the 10. 7 cm flux to the decay constants k for the upper, middle, and 
lower portions- of the topside ionosphere. Figure 9, therefore, shows the basis 
of obtaining the 3 independent slopes of the topside ionosphere as a function of 
f 0 F2, latitude, and solar flux. 

A further correlation to investigate the seasonal effects on k was carried 
out with some 15, 000 totally different Alouette soundings and fluctuations in the 
k values of i 15% were noted from the average spring and autumn values. The 
seasonal variation is monitored by observing the change in the daily maximum 
solar zenith angle from the equinoctial mid-day value. Figure 10 shows the 
seasonal fluctuation in k for each of the three layers in the topside profile. 
There is considerable evidence that this seasonal relationship has an added 
local time factor and this point will shortly be under investigation. 

Examination of the upper part of the’nose 1 of the .N-h profile is difficult 
because topside sounding information rarely gives any values in this region. 
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Evidence from many leading scientists also implies that the topside profiles 
have about a +4% error in the effective distance from the sounding satellite 
indicating the obtained topside profiles are too low near the peak. This 
evidence is based on comparisons with two-frequency data, backscatter 
results, Faraday rotation and overlap tests, etc. Preliminary results in 
this empirical model showed that a parabola in this region gave the better 
comparison with integrated total electron content when compared with two- 
frequency and Faraday rotation data. A simple parabola having a half thick- 
ness y t was fitted between the bi-parabola and the exponential layer. Upon 
initial test y^ was seu equal to the half . thi cknes s of the bi-parabola y n for £qF 2 
values below 10. 5 MHz, and y t increases with f Q F2 values rising above 
10. 5 MHz,. Further investigations of this problem are planned in future work. 

The final step in predicting the shape of the ionosphere is arranging for 
the gradient in the upper parabolic layer to be the same as the gradient in 
the lowest part of the topside exponential layer. This is the case at a distance 
d = _l/k [ (l+yf.k 2 )a -1] above the height of the maximum electron density. 

A, 1.3 Bottom side Ionosphere 

Modeling the bottoms} de ionospheric profile was a somewhat easier 
task because Tor each profile the value of f Q F2 was known and the electron 
density versus height profile from h nl to h ffiaJc was also known. Once more the 
geographic effect of longitude was eliminated and replaced with the more simple 
'local time correlation. From Figure 7 we see that the equation of the lower 
iayer is a parabola squared or a bi-parabola. This was found in general to 
fit the real profile somewhat better than a simple parabola. The unknown in 
this equation is the half thickness of the layer y a and in the reduction of the 
data the y B value was treated in a similar way to a topside k value. 

The irregularities in the ionosonde data due to the lower layers of the 
ionosphere were smoothed out because the prime objective of the work was to 
simplify the model, but keep the total content as accurate as possible. The 
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sounding data was therefore integrated up to the peak electron density (N B ) and 
forced to fit the bi -parabolic equation along with the value of N a obtained from 
the sounding. In each instance the value of y a was computed ready for further 
correlation. 

A number of real profiles from various stations at different local times 
were compared with the computed profile and excellent agreement found. 

A further 12, 000 soundings from all 14 stations were analyzed and the computed 
value of y m compared to the actual measured value. These results are shown • 
in Figure 11 along with the RMS errors. The two tests indicate that the bi- 
parabolic profile is, on average, in close agreement to the real profile. 
Investigations, similar to those carried out for the topside decay constants/ 
correlated y B with solar flux fo F2, local time and . season. Surprisingly, 
no direct correlation was found between y m and solar flux, but a definite 
correlation existed in local time and also in the solar zenith angle at local noon 
which represents the season. _ ' 

Figure 12 indicates how y ffl can be determined from local time and f Q F2, 
and Figure 13 shows the seasonal update as a function of local time for the 
sunrise, sunset, night and daytime period. In the cases where f Q F2 was 
larger than 10 MHz the local time curve fluctuated very little from the 10 MHz 
curve. All of the curves displayed have not been hand smoothed; due to the 
large data base the average of all values taken every hour fit precisely on 
the lines shown. 

The remaining unknowns which are needed to compute the profile are 
f 0 F2 and the height of that value; by far the most important of these being 
the value of f Q F2. 
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A. 1. 4 Predicting f 0 F2 


Severe horizontal gradients in 2 exist within the ionosphere as can 
be seen by examining Figure 14; In fact even if the value of * s known 

directly above a station, it can change considerably over the whole 'visible 1 
ionosphere from that site. Figure 14 is a predicted status of f Q F2 over the 
world at 6. 0 am during August .1968 and two types of severe gradients are 
immediately noticeable, one due. to sunrise causes rapid changes in f Q F2 m an 
east to west direction and the other situated around the equatorial anomaly 
occurs primarily during the afternoon and early evening and causes severe 
gradients in the north to south direction. Two hypothetical stations, A and B, 
are marked on Figure 14 along with the ionosphere ’visible' from those sites. 

In case A the value of f 0 F2 changes from 1 1.5 MHz directly overhead to 5 MHz 
on the southern horizon. This change must be squared when converting to 
electron content hence a difference of a factor of over 5 in the vertical content 
arises before correcting for elevation angle effects. Similar gradients exist 
over half the earth's surface at some time of the day and it is therefore 
imperative to model these gradients in any ionospheric model. 

For many years NOAA {formerly CRPL and 1TSA) have been engaged 
in the de velopment of numerical methods and computer programs for mapping 
and predicting characteristics of the ionosphere used in telecommunication s . 

The most advanced method for producing an f 0 F2 model undoubtedly comes 
from their work. Jones, Graham & Leftin (Reference 2 ) describe their 

techniques on how a monthly median of the F2 layer critical frequency (f 0 F2) 
was developed from an extremely large worldwide data base. In fact the gradient 
map shown in Figure 14 is a result of this work. We have already shown that it 
is important to include the horizontal gradients of fo^"^ an y analysis and the 
work by Jones et al is undoubtedly the only satisfactory approach to this problem. 

The document by Jones et al describing this work includes a Fortran 
program which, with monthly coefficients obtainable from NOAA, enables the 
monthly median value of £qF 2 to be computed above any point in the world at 
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any time. This program was primarily written to accept monthly coefficients 
using an average sunspot number, but more recent work by Jones & Obitts 
(Reference 3. ) has described a more generalized set of coefficients which 
provides annual continuity and uses more extensive analysis. These generalized 
coefficients can be obtained from the Ionospheric Prediction Services, NOAA, 
Boulder, for a sunspot number or a solar flux approach. The value of a monthly 
median f 0 F2 can be computed on a worldwide basis centralized around the specific 
day in question rather than the 15th of the month; it can also be based on a 
12-month running average of solar flux or sunspot number. Private communi- 
cation with Mrs. Leffcin at NOAA indicates that the solar flux approach is likely 
to provide more accurate values of f 0 F2 than the use of the sunspot number. 

For the ionospheric profile under discussion, it was decided to use the 
generalized f 0 F2 coefficients from NOAA incorporating solar flux thereby 
eliminating any need to purchase monthly data from them. The program was 
made self-contained and enabled a monthly median f F2 to be produced above 
any surface position for any time of day or season and any twelve month 
running average of solar flux. 

The question now arises as to how good these monthly median values 
are and how much error is introduced by day to day fluctuations. Many daily 
soundings were analyzed and the monthly median value computed; these were 
compared with the monthly median predicted values and the actual day to day 
fluctuations. Some typical results are shown in Figure 15. It is seen that the 
monthly median predicted values are indeed very close to the actual measured 
value, but the day to day fluctuations can be as large as +75%. A technique 
therefore had to be derived to bring the computed monthly . median value closer 
to the actual value. 

It would be pointless to use the daily value of solar flux in the generalized 
coefficient set which had been built up using a twelve month running average, 
but it was thought possible that there may be a relation between the difference 
in f 0 F2 from monthly median to daily value and the difference in the 12-month 
running average of solar flux to the daily value. 



Approximately 6, 000 real values of f D F2 from 13 stations widely spread 
in latitude, longitude, and solar cycle were compared with the predicted values 
using the NOAA solar flux method. A very surprising result emerged and can 
be explained by referring to Figure l6. Eliminating the data from stations 
close to the magnetic poles which did not quite follow the trend of the other 
stations a comparison between the difference in daily and 12-month flux value 
and the percentage difference of computed and measured f 0 F2 showed all stations 
having a very similar bias. Figure 16. shows this comparison where the stations 
having similar latitude were averaged quoting their mean magnetic latitude. The 
fact that the lines did not pass through the zero points in the graph undoubtedly 
indicates an erroneous bias in the NOAA predictions, but results help one to 
update substantially the monthly median f F 2 value on a daily basis. Further 
comparisons were carried out with two years of hourly f 0 F2 values obtained near 
solar maximum from Hawaii and the results fit perfectly in the latitude position 
expected in Figure 16. By these means it is possible to come somewhat nearer 
ine actual aaiiy Value of f Q F2. Further accuracy can be. derived by update 
from stations within the general area if this is available and the investigation 
of this approach will now be explained. 

In order to investigate the size of an area from which ionosphe ric value s 
would show similar deviations from normal, many comparisons of three or 
more stations were investigated for random dates. It is well known that 
magnetic disturbances • can effect the ionosphere above one station in one 
direction and a nearby station in an opposite direction. For this reason 
investigations of disturbances were not carried out near to the magnetic poles. 
Over 100 groups of stations from various continents and having similar 
longitudes were compared in similar ways. Figure 1 7 is a typical result of 
such a test and shows f Q F2 disturbances being recorded simultaneously at 
sites 1, 000km apart. The percentage error in the predicted f Q F2 value when 
compared to the real value was noted to be similar in 90% of the cases where 
stations were within 2, 000km of one another in a longitudinal direction and 
investigations over the ’quiet’ North American continent show improvement 
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in 9 out of 10 cases when f 0 F2 was updated with information from across the 
continent; or 3, 000 to 4, 000km. However, in general, the update procedure 
is restricted to information from within 2, 000km of the evaluating station. 

A i. 5 Predi cting the Height of the Maximum Layer 

' In order to predict the real height of f 0 F2 the M(3000)F2 predictions from 
NOAA were used. To explain the terminology: 

M(3000)F2 = M FACTOR = MUF(3000)F2 / f 0 F2, 

where MUF(3000)F2 is the maximum usable frequency to propagate by 
reflection from the F2 .layer a distance of 3, 000km. The M(3 000)F2 predictions 
can be calculated on a monthly basis from a generalized set is sued by NOAA 
and provide the monthly median value as a function of sunspot numDer. 

Knowledge of this factor along with the f Q F2 value enables the height of 
the layer to be calculated using the equations of Appleton & Beymon (Reference 1 
If M is the M(3 000)F2 factor and one assumes that y m divided by the height of 
the bottom edge of the lower layer is greater than 0. 4, then it is possible to 
derive the following polynomial, 

h = 1346. 92 - 526. 40M + 59. 825M 2 ., 

B 

where h^ is the required height.. 

A. 2 Model Accuracy : • 

As a means of testing the accuracy of the model, an intense comparison 
with Faraday rotation data has been performed as well as tests with two 
frequency data, actual ionospheric profiles, and use in orbit determination 
programs. 

Remarkable improvements have been noticed in precise orbit determination 
systems and the model has reduced the number of iterations needed for the 
program to converge as well as the size of the residuals by up to a factor of 
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four. Excellent results have been noted with orbit programs using elevation 
angle, range and range rate systems. 

The most extensive tests were carried out by comparing Faraday rotation 
data for seven stations from Hawaii to Puerto Rico to ALaska looking at the 
ATS1, ATS3, and SYNCOM3 satellites. In all, over 100 station months of 
continuous data were used during the years 1965 and 1967-1969 with data 
taken every hour. The integrated model data was compared with these actual 
results; update situations were also investigated. The results are shown in 
Figure 18 where the percentage of the ionosphere removed with the model 
is shown. In general, between 75 and 90% of the ionospheric effects are 
removed and these circumstances are for solar maximum conditions. 
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APPENDIX B 


B. 1 Earth's Magnetic Field Model 


The model computes the earth's magnetic field components at a 
desired location following the spherical harmonic analysis of the 
magnetic field by Chapman and Bartels (Reference 1) and using the 
coefficients g“ , h“ given by Jensen and Cain (Reference 2) for Epoch I960. 
The X-north, Y-east, and Z- vertical (up) components of the magnetic 
field are computed for any location, defined by its latitude 0, longitude X , 
and height h above the earth's surface. Introducing the colatitude 
cp=90° -0 and the ratio R=R e /(R e +h), where is the radius of the earth, 
the components X, Y, Z are given by. 


6 n — “I'* 

x =I { r " + M jL b (cosep) [^g® cos(mX) + h° sin(mX ) Jj 

— 1 a 

Y = J_ V ( m „ (coscp) [ g ; sin(rra) - h; cos(mX) J) 

Z.T |(n+l) £ P n , H (cosep) [g* cos(mX ) + h® sin(m\ ) J } 

b=u b=o 

The multiple of the associated Legendre function is given by, 


„ _ (n-m )(n-m- 1 ) n -w-p 

Pfl» n ( COSC P) ~ sin® CP L COS CP - cos $ 


i n CD L ' 

• L 2(2n_l) 

(n-m)(n-m- I )(n-m-Z )(n-m-3 ) cos b- d -4 1 

(2) (4) ( 2n- 1 ) (2n- 3 ) <P ' J 
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